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ABSTRACT 



A particle acceleration mechanism by radiation pressure of precursor waves in a 
relativistic shock is studied. For a relativistic, perpendicular shock with the upstream 
bulk Lorentz factor of 71 S> 1, large amplitude electromagnetic (light) waves are known 
^' to be excited in the shock front due to the synchrotron maser instability, and those waves 

I ■ can propagate towards upstream as precursor waves. We find that non-thermal, high 

2 ■ energy electrons and ions can be quickly produced by an action of electrostatic wakefields 

, generated by the ponderomotive force of the precursor waves. The particles can be 

quickly accelerated up to emax/Ti^T-eC^ ~ 71 in the upstream coherent wakefield region, 
and they can be further accelerated during the nonlinear stage of the wakefield evolution. 
\ The maximum attainable energy is estimated by Smax/Ti^ieC^ ~ Lsys/{c/ujpe), where 

Lsys and c/ujpe are the size of an astrophysical object and the electron inertial length. 



' respectively. 



Subject headings: acceleration of particles — radiation mechanisms: nonthermal 



, cosmic rays — plasmas — shock waves 

> 

X ■ 1. INTRODUCTION 

^' 

Acceleration of particles is a universal phenomenon occurring in a variety of objects such 
as solar flares, interplanetary shocks, the heliospheric termination shock, stellar flares in young 
stars, astrophysical jets, pulsars and black holes (e.g., Begelman et al. 1984; Asseo and Sol 1987; 
Meisenheimer et al. 1989; Koyama et al. 1995: Carilli and Barthel 1996; Junor et al. 1999). After 
the discovery of cosmic rays in the early last century by the balloon experiments by Hess (1912), 
the origin of the cosmic rays has been a long standing problem in astrophysics and cosmic ray 
physics. Particle energies up to 10^^'^ eV are thought to be generated by supernova shocks, and 
the energies of 10^^'^ eV and beyond have been suggested to be of extra-galactic origin such as 
active galactic nuclei (AGN) jets and gamma ray bursts (GRBs) at cosmological distances (e.g., 
Thompson 1994; Vietri 1995). The black holes and pulsars may be candidates as well. Yet it is not 
well understood how these powerful astrophysical objects can produce such high energy particles 
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(e.g., Hillas 1984; Cesarsky 1992). One of the goals of particle acceleration is to explain the origin 
of those high energy cosmic rays in the universe. 

The generation processes of non-thermal particles may be categorized into three types: One 
is the energy conversion from the magnetic field towards non-thermal particles, the second type is 
the energy conversion from the bulk flow energy, and the third one is the energy conversion due to 
large amplitude waves. Magnetic reconnection is one of the examples of the magnetic field energy 
dissipation processes, and particle acceleration in magnetic reconnection was proposed to operate 
in various astrophysical sites such as radio jets, pulsar winds, and soft gamma repeaters (e.g., 
Coroniti 1990; Romanova and Lovelace 1992; Thompson and Duncan 1995; Kirk and Skjaeraasen 
2003). Recently several people discussed that a non-thermal power law energy spectrum can be 
generated by relativistic reconnection, where the Alfven speed is close to the speed of light (e.g., 
Zenitani and Hoshino 2001, 2005; Jaroschek et al. 2004). The importance of reconnection for 
astrophysics is that it may occur wherever there is a shear of magnetic field. 

The second type of the non-thermal particle acceleration is the energy conversion from low- 
entropy bulk flow energy to non-thermal particles. A shock wave is known to be one of the most 
efficient energy conversion processes in this class. Wherever the supersonic flow exists, the genera- 
tion of the shock wave is an inevitable consequence when the flow encounters with a surrounding 
plasma medium as an obstacle. Most of the bulk flow energy can be converted into not only thermal 
energy but also non-thermal particle energy. One of the accepted schemes of the shock acceleration 
is the so-called diffusive shock acceleration/Fermi acceleration (e.g., Axford et al. 1977; Bell 1977; 
Blandford and Ostriker 1978; Drury 1983), in which the particle is accelerated by repeated cross- 
ing of the shock front between upstream and downstream turbulent states. The most important 
property of diff'usive shock acceleration (DSA) is a universal power-law energy spectrum with the 
power-law index of 2, and the non-thermal energy spectrum can be determined by the density com- 
pression of plasma between the shock upstream and downstream, which does not strongly depend 
on the detailed plasma parameters. By virtue of this behavior, DSA has been widely accepted 
as the most plausible shock acceleration mechanism for the application to astrophysical situations 
(e.g., Blandford and Eichler 1987; Jones and Ellison 1991). 

In addition to DSA, several acceleration schemes with fast acceleration mechanisms have been 
discussed in the context of direct acceleration, in which the particle can be accelerated continuously 
by traveling along the electric fleld direction. The shock drift acceleration is one of the schemes 
of direct acceleration, which operates at the shock front (e.g., Sonnerup 1969; Armstrong et al. 
1985). A particle can gain its energy during the so-called gradient B drift motion parallel to the 
motional electric field at the shock front. The shock surfing acceleration (SSA) is also another 
direct acceleration scheme by using the process of particle trapping by the electrostatic/magnetic 
potential well at the shock front (e.g., Sagdeev and Shapiro 1973; Katsouleas and Dawson 1985; 
Hoshino 2001). SSA is recognized as highly efficient non-thermal particle acceleration mechanism, 
and is believed to serve at least as "seed" particles which are subject to further acceleration to 
much higher energies by DSA (e.g., Dieckmann et al. 2000; McClements et al. 2001; Hoshino and 
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Shimada 2002; Amano and Hoshino, 2006). 

The third type of particle acceleration is to utilize the wave energy as free energy source (e.g., 
Robinson 1997). MHD turbulence is often discussed as an elementary process of particle acceleration 
in various situations (e.g., Bieber ct al. 1994), and the resonant scattering of particles with the 
turbulent MHD waves is known to contribute to the non-thermal particle acceleration as well as 
the plasma heating. (Since the particle acceleration in DSA comes from the scattering of particles 
within the MHD turbulence, our classification is just for the sake of convenience.) In addition to 
the turbulent acceleration, a large-amplitude coherent wave can be invoked for non-thermal particle 
acceleration, because a wave packet of a large-amplitude wave exerts pressure on plasma through 
the so-called ponderomotive force (e.g., Kruer 1988; Mima and Nishikawa 1988; Esarey et al. 1996; 
see also Appendix A). In this paper, we argue a relatively new acceleration process by the radiation 
pressure of large-amplitude electromagnetic waves excited in a relativistic shock wave. We focus 
on a large-amplitude electromagnetic wave generated by the collective plasma process, and discuss 
how the non-thermal particles can be effectively produced during the dissipation of the wave energy. 

The above idea that large-amplitude electromagnetic waves in a relativistic shock are responsi- 
ble for high-energy particle acceleration is not necessarily new. Recently Lyubarsky (2006) discussed 
an energy conversion process from ion to electron throught the upstream electromagnetic precur- 
sor wave in a relativistic shock wave. It is known that the large amplitude precursor waves can 
propagate towards upstream from the shock front, and he demonstrated that the velocity of the 
electron guiding center decreases inside the large amplitude waves, and that the electrostatic field 
generated by the relative motion between ions and electrons acts on the particle acceleration. 

We extend the intriguing study done by Lyubarsky (2006) and demonstrate that a wakefield 
acceleration mechanism can operate in the relativistic shock upstream region. After Tajima and 
Dawson (1979) first proposed the idea of the wakefield acceleration in an ultra- intense laser pulse, 
the particle acceleration by the wakefield has been widely investigated in laboratory laser plasma 
experiments (e.g., Easrey et al. 1996; Mouron et al. 2006). It is discussed that some electrons can be 
preferentially accelerated by the electrostatic wakefield generated by the intense laser pulse. Based 
on this idea, the wakefield acceleration of Alfvenic waves in astroplasma settings has been argued 
by Chen et al. (2002), who suggested that ultra-high energy cosmic rays may be generated by the 
wakefield acceleration. We investigate an intrinsic wakefield acceleration in relativistic shocks, and 
show that the electromagnetic (light) precursor wave self-consistently generated in the shock plays 
a unique role on particle acceleration, which has some diff'erent aspects from the laboratory laser 
experiments. By the action of the ponderomotive force to electrons due to the large-amplitude 
electromagnetic precursor wave, electrons can be expelled from the precursor wave region, and 
then the large amplitude electrostatic wakefield behind the packet of the precursor wave can be 
generated, which in turn can accelerate some electrons and ions towards ultra-relativistic energies. 
By using a particle-in-cell (PIC) simulation, we demonstrate that the wakefield acceleration provides 
a strong electron and ion acceleration in front of the relativistic shock waves. 
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2. 



RELATIVISTIC PERPENDICULAR SHOCK STRUCTURE 



2.1. 



Review of Relativistic Shock Simulation 



The dynamic structure of the relativistic perpendicular shock and its particle acceleration has 
been studied by many people by means of the particle-in-cell (PIC) simulations (e.g., Langdon et 
al., 1988; Gallant et al. 1992; Hoshino 2001; Nishikawa et al. 2003; Hededal et al. 2004). It 
was discussed that strong plasma thermalization occurs in a relativistic shock front region, where 
a rich repository of plasma instabilities is taking place. The basic structure of the relativistic 
perpendicular shock in a pair plasma is controlled by the so-called a± parameter, which is defined 
by the ratio of the upstream Poynting flux to the upstream particle-born energy flux, namely, 
(T-t = Bf /{STrNi^imeC^), where Bi, Ni and 71 are the upstream magnetic field, the upstream 
density, and the Lorentz factor, respectively. The a± parameter is Lorentz invariant. 

For the case of c-i- < 1, it is found that the plasma heating quickly occurs through electromag- 
netic and electrostatic waves generated at the shock front. The origin of those waves is understood 
by the synchrotron maser instability of the gyrating particles at the shock front, which is pro- 
duced by the interaction between the shock front and the incoming particles from the upstream 
region (Hoshino and Arons 1991). Those waves generated at the shock front can propagate into 
both downstream and upstream, and the upstream propagating electromagnetic waves form the 
precursor waves. The upstream precursor waves are a unique feature of the relativistic shock. It is 
discussed that the amplitude of the precursor waves normalized by the upstream ambient magnetic 
field increases with decreasing a± , but the emission of the precursor wave does not affect the shock 
dynamics for a relativistic shock in electron-positron plasmas. 

After the quick thermalization of the upstream, low-entropy plasma at the shock front, a 
relativistic hot plasma is formed in downstream. The energy spectra in downstream can be approx- 
imated by a relativistic thermal-Maxwellian, and did not show any evidence of non-thermal particle 
acceleration for a shock with pure electron-positron plasmas (e.g., Langdon et al. 1988; Gallant 
et al. 1992), while Hoshino et al. (1992) demonstrated that significant nonthermal electrons and 
positrons can be produced in the shock downstream, if ions with a larger inertia contaminated in 
the electron-positron plasma. They discussed that the gyrating ion energy near the shock front 
can be quickly transferred into the positron and electron non-thermal energy through the relativis- 
tic cyclotron resonance with the higher harmonic waves generated by the ion synchrotron maser 
instability. 



In this paper, contrary to the previous shock studies of the pure pair plasma or of the pair 
plasma with a minority ion population, we discuss the ion-electron shock. We show that the 
precursor waves play a very important role for the electron and ion acceleration. We use a one- 
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Ion-Electron Shock Simulation 
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dimensional, particle-in-cell (PIC) simulation. In the simulation study, the geometry of the shock 
is the same as for the previous shock simulations in pair plasmas. A low-entropy, relativistic 
plasma consisting of electrons and ions is injected from the left-hand boundary region which travels 
towards positive x. At the particle injection boundary, the plasma carries a uniform magnetic 
field Bz and a uniform motional electric field Ey = x B^/c, polarized transverse to the flow Vx- 
The downstream right-hand boundary condition at x = is a wall where particles and waves are 
reflected. Under the interaction between the plasma traveling towards positive x— direction and 
the reflected particles and waves from the right-hand boundary, the shock wave is produced, and it 
propagates backward in the —x direction. After the formation of the shock front, the self-consistent 
energy dissipation is provided through a collective plasma process at the shock front. The grid size 
is set to be O.OSSOc/wpe) and includes 100 particles/cell for each species at the left-hand injection 
boundary. The plasma frequency LOpe = \/ AttNi e^/jiuie = \/ ^imie^ / rUe, where Ni and ni are the 
plasma density in the simulation (shock downstream) frame and in the proper frame (upstream 
frame), respectively. The speed of light c is normalized to be unity. The other plasma parameters 
are as follows: the upstream Lorcntz factor for the bulk flow speed 71 = 10 and the mass ratio of 
ion to electron mj/me = 50 in order to save computational CPU time. The a parameter is set to 
be (Ttotai = 5iV(47rAri7i(mj + m^y) = 2 x 10-^ and = Bf/{4TrNijimec'^) c± 0.1. 

Shown in Figure 1 is a snapshot of the relativistic, perpendicular shock at tcope = 5431 obtained 
in our PIC simulation. The x-axis is normalized by the relativistic electron inertial wavelength, 
2'rrc/LOpe. Prom the top, the ion phase space Uix in the x direction, the ion phase space Uiy in 
the y direction, the electron phase spaces U^x and U^y, the longitudinal electric field Ex, and 
the transverse magnetic field B^. These physical quantities are normalized by the corresponding 
upstream parameters. In the phase space plots, the color contour means the phase space density (the 
reddish/bluish region is high/low phase space density). A well-developed shock front can be formed 
at the center of the simulation box, i.e., x/{2Trc/L0pe) ~ —400. The right-hand region corresponds 
to the shock downstream occupied by relativistically hot plasmas, while the left-hand region is 
the shock upstream/transition region with super-sonic flows. We can see that large amplitude 
electromagnetic waves (Bz) generated at the shock front are propagating into the shock upstream 
region. The precursor wave B^ is associated with the transverse electric field Ey (not shown here). 
At the time of tupe = 5431, the tip of the precursor waves reaches x/{2Trc/ujpe) ~ —850, and 
we find that the propagation velocity of the precursor wave is about (850 x 27rc)/5431 ~ 0.98c, 
which is almost equal to the speed of light. In contrast with non-relativistic shock waves, the 
relativistic shocks can transfer the upstream kinetic energy not only into lower frequency MHD 
waves but also the higher frequency electromagnetic waves. We will discuss later the high frequency 
electromagnetic wave energy density contains about 10% — 100% of the upstream electron bulk flow 
energy density. 

In addition to the large amplitude precursor waves of the transverse magnetic field B^, the 
longitudinal electric field Ex and the pre-heating of electrons can be clearly seen in the shock 
transition region. The longitudinal electric field Ex has much longer wavelength than the transverse 
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electromagnetic field B^. The modulation of the electron phase space Ugx is associated with the 

oscillatory electric field Ex- Within the oscillatory motion electrons are accelerated up to Uex/Ui < 
20. The key point of this paper is the pre-heating and pre-acceleration of electrons in the shock 
transition region. As we will explain later, the longitudinal electrostatic waves Ex can be understood 
to be the so-called "wakefields" excited by the high-frequency precursor electromagnetic waves B^. 

Let us carefully study the wakefield structure. The modulation of the electron phase space 
in Uex has a good correlation with the wave phase of Ex at the front side of the wakefield. The 
negative region of Ex corresponds to the electron accelerated region. As approaching towards the 
shock front region, the correlation between the electron phase space and the wakefield is smeared 
out, and they show a complex structure. It is interesting to note that many electrons are accelerated 
up to several tens of the upstream bulk flow energy, namely the maximum energy of the accelerated 
electrons is larger than the upstream bulk flow energy of ions, because in our simulation the mass 
ratio of ion to electron is assumed to be 50. In addition to the electron acceleration, we can see 
also the acceleration/heating of ions inside the upstream wakefield region. 

Shown in Figure 2 are the ion and electron energy spectra in the shock upstream and down- 
stream. The upstream and downstream spectra are denoted by the dashed lines and the solid 
lines, respectively. Both upstream and downstream spectra are obtained by integrating the particle 
energy measured in the shock downstream (laboratory) frame. Some electrons in upstream are 
already accelerated up to e/71 mgC^ ~ 215, and the maximum energy of electrons in downstream is 
almost same as the upstream maximum energy. The flux of high energy electrons in downstream 
is enhanced due to some plasma acceleration and heating processes in and around the shock front. 
The acceleration of ions is not significant compared with electrons, but the maximum ion energy 
of jruiC^ with 7/71 ~ 8 is of the same order of the electron maximum energy, because the mass 
ratio of ion to electron used in our simulation is set to be 50. The ions in upstream show a narrow 
spectrum compared with the electrons, but we find some ions accelerated up to 7/71 ^ 6. By 
taking into account a slow response of ion acceleration, we think that both ions and electrons are 
equally accelerated by the electrostatic wakefields. 



3. WAKEFIELD GENERATION BY PRECURSOR WAVES 

In the PIC simulation, we found that the upstream particles are accelerated through the 
generation of the electrostatic field Ex- As shown in Figure 1, the large-amplitude, precursor 
electromagnetic wave Bz in upstream seems to trigger the electrostatic field Ex- Let us discuss how 
and why the electrostatic field Ex can be generated in the shock upstream. 



- 7- 



3.1. Electrostatic Wakefield in Relativistic Shocks 

After the first discovery of tiie wakefield acceleration by Tajima and Dawson (1987), great 
progress of the wakefield acceleration has been made in theory and laboratory laser experiments 
(e.g., Sprangle et al. 1990; Esarey et al. 1996). The conversion of the energy of electromag- 
netic waves (photons) into intense electron beams has been extensively discussed in laser-plasma 
interactions, as an important application to laser-plasma accelerators and the fast igniter fusion 
concept. 

If a packet of large amplitude electromagnetic waves is injected into the plasma medium, 
the electromagnetic wave can expel electrons in the front of the wave packet so as to induce the 
longitudinal polarization electric field in order to maintain charge neutrality. Then the polarization 
electric field acts on the longitudinal electron motion, which in turn leads to plasma oscillations. 
In this way, an electrostatic Langmuir wave can be excited behind the electromagnetic wave packet 
(e.g., Kruer 1988). 

This interaction of the electromagnetic wave with the Langmuir wave can be understood by 
the so-called ponderomotive force (cf. Appendix A), and the force -Fpond is written in the shock 
upstream frame where the upstream plasma flow is at rest, 

-^pond = -eV^pond, (1) 

where ^ 

0pond = 7 — r2-|-SoP- (2) 
4 mo^Q 

Eq and Uo are the amplitude of the injected electromagnetic (light) wave and its frequency, re- 
spectively. The physical quantities in the shock upstream frame are denoted with the symbol 
of "tilde", while those in the laboratory frame where the downstream plasma flow is at rest are 
denoted without "tilde" . The ponderomotive force is proportional to the gradient of the wave pres- 
sure, and is independent of the sign of the charge. However, the ponderomotive force for electrons 
is much stronger than that for ions. By neglecting the ion response to the ponderomotive force, 
the induced electrostatic field -E^ake can be estimated by balancing the electrostatic force and the 
ponderomotive force, 

-^wake ~ ~-Epond- (3) 

As we will discuss later, the amplitude of the electromagnetic precursor wave becomes large 
for the case of a relativistic shock. In this case we have to use the generalized ponderomotive force 
(e.g., Bauer et al. 1995), 

^pond = mc^^l + rjal, (4) 

where the normalized amplitude of the electromagnetic wave is given by do = eEo/meOOo = 
eAo/mcP and r/ represents the wave polarization: ry = 1 for circular and 1/2 for linear polar- 
ization. The normalized amplitude oq is Lorentz invariant. For the case of Oq = ao <C 1, we can 
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easily confirm tliat the generalized potential of Eq.(j3|) is consistent with the potential in a small 
amplitude regime of Eq. . 

3.2. Scale Length of Wakefield: Raman Scattering and Self-Modulation 

The ponderomotive force works onto the plasma under the gradient of the wave amplitude. As 
shown in Figure 1, the ponderomotive force is effective not only at the tip of the precursor wave 
in upstream but also inside the precursor wave, because the wave amplitude inside the precursor 
region can be periodically modulated. We think that one of the modulation processes comes from 
the so-called Raman scattering instability in the collective plasma process, and the other is the self- 
modulation due to the nonlinear interaction between the accelerated electrons and the precursor 
wave. 

Let us first quickly review the stimulated Raman scattering process. As we mentioned the basic 
mechanism of the Raman scattering in Appendix B, a large-amplitude electromagnetic (light) wave 
can decay into two waves: one is a longitudinal Langmuir wave, and the other is a scattered 
electromagnetic (light) wave. Let us assume the electromagnetic wave is propagating into a plasma 
with a density perturbation 6n. Since the phase speed of the light wave becomes larger in the denser 
region of the density perturbation, the refractive index becomes large and the light wave can be 
pushed out of the positive density perturbation region into the negative density perturbation one. 
The non-uniform light wave generates the ponderomotive force, which in turn enhances the initial 
density perturbation, which is organized by the Langmuir wave. According to Raman scattering 
theory, the wavelength of the Langmuir wave is scaled by the electron inertial length c/wpe in the 
plasma frame, i.e., the shock upstream frame, if the frequency of the electromagnetic wave is much 
larger than the plasma frequency ujpe- 

The periodic variation of the precursor wave amplitude can be also initiated by a self-modulation 
effect between the electrons accelerated by the precursor wave and the precursor wave generated by 
the electron flow. Consider the early phase of the shock evolution when the precursor wave is just 
launched from the shock front. In the leading edge of the precursor wave, the wave amplitude is 
not necessarily large, but a finite electric field can be generated by the ponderomotive force of 
the precursor wave. The polarization electric field is negative in the tip region of the precursor 
wave, and an electron can be accelerated towards the positive x direction. The accelerated electron 
now interacts with the shock front, which in turn can emit a larger amplitude precursor wave by 
the synchrotron maser instability at the shock front. Therefore, the polarization electric fields can 
be enhanced by this feedback effect. During the feedback process, electrons are accelerated against 
the induced electric field, and the plasma oscillation is initiated. The above feedback process will 
couple to the Raman scattering, and leads to a self-modulation of the precursor wave. 

Let us discuss the scale length of the wakefield perturbation observed in our simulation. Since 
the scale length of the wakefield L is observed in the laboratory/simulation (or the shock down- 
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stream) frame, we have to take into account the Doppler shift effect on the wakefield in order to 
compare to the above Raman scattering theory. Since the wakefield is propagating towards the 
upstream direction in the shock upstream frame, the wavelength of the wakefield in the labora- 
tory frame should be elongated. By virtue of the Lorentz transformation, the frequency uj and its 
wavenumber k in the shock downstream (laboratory) frame can be related to io and k in the shock 
upstream frame, 

/ a;/c \ ^ / 71 ^^1 ] [ ^/c \ 
\ k I )\ k ) 

where ui = y^^f — 1 = jiPi- For the case of the scattering in our simulation, we have u) ~ ojpe and 
k ~ —ujpe/c. Therefore, the wave number k in the laboratory frame is given by, 

fc=.7i(/3i-l)^--^^, (5) 
c c 271 

and the scale length of the wakefield structure L in the laboratory frame can be evaluated by 

L~27i — . (6) 

The above estimate can be compared with our simulation result. Figure 3 shows the wave 
power spectra in two different regions of the relativistic shock upstream at tiVpe = 5432 when the 
shock wave has been well developed. Regions (a) and (c) are in the tip of the precursor wave, while 
regions (b) and (d) are inside the precursor wave train. The wave number is normalized by the 
wave number of the electron inertial length kp = uJpe/c. The electrostatic wave spectrum has the 
maximum peak amplitude around k/kp ~ 0.05 at the tip of the precursor wave. Since we used the 
upstream bulk Lorentz fact of 71 = 10 in the simulation, the observed wave number k/kp ^ 0.05 is 
expected to be equal to (271)^^ = 0.05 from Eq.([5]). We find a very good agreement between the 
simulation and the theory of the Raman scattering process. 

Inside the precursor wave region, those waves show broad band spectra due to some nonlinear 
wave instabilities, and the peak amplitudes of both electromagnetic wave and the wakefield are 
slightly shifted towards longer wavelength regions with k/kp ~ 2 and k/kp ~ 0.03, respectively. 
As we will discuss later, the longer wavelength of the wakefields can slightly enhance the particle 
acceleration efficiency. 



3.3. Amplitude of Wakefield 

Since we obtain that the scale length of the wakefield is about the electron inertial length c/ujpe, 
the amplitude of the electrostatic wakefield i?wake discussed in Eq.([3]), which is Lorentz invariant, 
can be given by using the quantities in the downstream frame, 

TjO?' 

E-wake — -^wake — / ^ ^wb ; C^) 
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where -E^b = ^eC^pe/^ is the so-called wave-breaking field, in which plasmas are capable of sup- 
porting large amplitude electrostatic waves with phase velocities near the speed of light. We have 
neglected factors on the order of unity. We may rewrite the amplitude of the electrostatic wakefield 
normalized by the upstream magnetic field as, 



-E'wakc / 1 



(8) 



71, 



where we have assumed the electron cje = / {Att Ni^irrieC^) and Bi ~ Ei. In order to estimate 
the amplitude of the wakefield in Eq.dS]), we have to know further the normalized amplitude ao = 
eEo/meCijQ, namely, the amplitude of the precursor electric field £"0 and its wave frequency ujq. 

The dynamic behavior of the precursor wave has been discussed by Gallant et al. (1991). They 
investigated the structure of the relativistic shock waves in pair plasmas by using a one-dimensional 
PIC simulation code, and discussed in detail about the emission process of the precursor waves as a 
function of a± = i?^ /(SvrA'^i 71 mgC^). They showed that the frequency of the electromagnetic wave 
loq in the downstream frame can be given by 

— ~ Tshock ( Vo-± + 2 + /3shock Vo-± + 1) , (9) 



where 7shock = l/y 1 ~ /^shock Lorentz factor of the shock propagation speed. From the shock 

jump conditions, i.e., Rankine-Hugoniot relations, we obtain /3shock = T — 1 for u-i- ^ 1, where T is 
the ratio of specific heat. Therefore, we have 

Wo 



IjJpe 



2-3. (10) 



Their simulation study is performed for a pair plasma shock, but the above study is applicable 
for the ion-electron shock, because the emission process of the precursor wave is mainly controlled 
by an anisotropic velocity distribution function of electron/positron. The upstream cold incoming 
ions and electrons start a gyro-motion at the shock front, and their distribution functions become 
an anisotropic ring-type distribution perpendicular to the ambient magnetic field. Both ions and 
electrons are capable to emit collective electromagnetic waves through the synchrotron maser insta- 
bility (Hoshino and Arons, 1991). The electron ring-type distribution can generate a high-frequency 
X-mode wave, whoes phase velocity is almost equal to the speed of light, and the excited X-mode 
wave can propagate towards upstream in the form of the precursor wave. The ion ring-type dis- 
tribution, however, mainly generates the magnetosonic wave whoes phase speed is slower than the 
shock propagation speed, and the magnetosonic wave cannot propagate towards upstream. There- 
fore, the precursor wave in the shock upstream is controlled by the electron physics. In fact, as seen 
in Figure 3, the peak amplitude of the electromagnetic wave spectrum is k/kp ~ 3. By assuming 
that oj ^ kc for k ^ kp, we find that the simulation result agrees very well with the theoretical 
estimation of Eq. (fTOl) . 

Next we discuss about the amplitude of the precursor electromagnetic wave. The free energy 
of the generation of the precursor wave is the upstream bulk flow energy, and the wave energy 
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density would not exceed the bulk flow energy density. It is not easy to estimate analytically the 
energy conversion rate from the bulk flow energy to the precursor wave energy in our knowledge of 
plasma physics. Therefore we refer the estimation to the numerical simulation experiment. Figure 
4 shows the precursor wave energy Bq normalized by the upstream ambient magnetic field energy 
and by the upstream bulk flow energy, which is reproduced from the simulation result by Gallant et 
al. (1992). The precursor wave power normalized by the upstream magnetic field energy increases 
with decreasing the a±, while if it is normalized by the upstream bulk flow energy, we find that 
the wave energy density is almost constant with 10% for the range of 10~^ < cr± < 1. Namely the 
energy conversion rate from the upstream bulk fiow energy into the precursor wave energy econv is 
given by 

econv = B^/{8TrNi^imeC^) ~ 0.1. (11) 

For a large cr± > 1, the error bar becomes large, because the shock Mach number decreases with 
increasing a±, and because the shock wave is not necessarily well-developed yet. For a small 
(Te < lO"'^, not only the synchrotron maser instability but also the Weibel instability may play 
an important role on the shock dissipation process (e.g., Yang et al. 1994; Nishikawa et al. 2003; 
Hededal et al. 2004; Kato 2005). The emission of the precursor wave and the generation of the 
Wakefield should be carefully studied for such a small cJe regime. 

Based on the above discussion, the amplitude of the precursor electromagnetic wave in the 
ion-electron shock can be given by, 

En Bo /2enonv\^^^ . ^ 



Bi Bi V 

where we have used the electron cJe, which differs from a± by the factor of 2, and the normalized 
precursor amplitude oq is written in the downstream frame as, 

UJpe Eq ^pe j / \ 

OO - IW/Oe — =^ 71 v2econv ~ TlV^conv (13) 

One can find that the normalized amplitude becomes much larger than unity for a super-relativistic 
shock wave with 71 » 1, and the wakefield Ermwake in Eq.® becomes -Ewakc ~ ^0 for 7i ^ 1- 



3.4. Comparison of Amplitudes Between Theory and Simulation 

We compare the above theoretical estimations with our simulation results. In addition to 
(o"totah f''e) = (2 X 10^'^, 0.1) presented in Figure 1, we performed other simulations with the param- 
eters of (o"totahfe) = (8 X 10~^,0.4) and (3.2 x 10~^, 1.6). Except for a, we used the same plasma 
parameters as Figure 1, and we assumed 71 = 10 and the mass ratio of ion to electron is 50. The 
grid size A normalized by the relativistic electron inertial length c/ujpe is set to be 0.0389. The 
shock structures of those higher a shocks and their nonlinear dynamic evolutions were basically 
same as that in Figure 1, namely, the large amplitude, electromagnetic precursor waves emanate 
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from the shock front, and the electrostatic wakefields are generated behind the precursor waves. 
The strong modulations of the incoming electrons are observed inside the wakefield regions. The 
small deformations of the incoming ions are also seen as well. 

Shown in Figures 5a and 5b are the comparison of theory and simulation. Figure 5a shows 
the amplitudes of the electromagnetic precursor waves Bq as the function of the electron ag , while 
Figure 5b is the corresponding electrostatic wakefields -Ewake- In Figure 5a, the closed squares are 
the simulation results measured in the leading edge of the precursor wave, while the closed circles 
denote the maximum amplitude of the precursor wave obtained inside the precursor wave. Two 
solid curves are the theoretical estimation based on our discussions in Eq. (jl2p . The lower curve 
is estimated by assuming £conv — 0.1, while the upper curve is Cconv 

= 1.0. The amplitude of the 

precursor wave in the leading edge is almost consistent with the theoretical curve of econv = 0.1, 
while the maximum amplitude inside the precursor waves can be roughly scaled by econv = 1-0. 

As we have already discussed for the self-modulation of the precursor waves, the feedback 
process works between the modulated electron flow by the precursor wave and the enhancement 
of the precursor wave due to the accelerated electron flow. This feedback process can enhance the 
energy conversion rate Cconv) and it may saturate when the excited wave energy density become of 
the order of the bulk flow energy density. From the above discussion, we may adopt that 

econv ~0(1), (14) 

inside the precursor wave region in the downstream frame. 

In Figure 5b, the closed circles are the simulation results of the electrostatic wakefields -Ewakei 
which are found in association with the large amplitude electromagnetic precursor wave. The solid 
line stands for the theoretical estimation in Eq. ([8]). We have assumed that e^nv = 1-0, because the 
maximum amplitude of the precursor waves are modeled by e^nv = 1.0 from the above discussion. 
We have also used ao = 71 V 2econv ("^pe/ 1^0 ) with LOo/u)pe = 2.5. We find that the agreement between 
the theory and the simulation is very good. 

4. PARTICLE ACCELERATION BY WAKEFIELDS 

Let us discuss electron acceleration by the stimulated wakefield in the region of the upstream 
precursor wave. As seen in Figure 1, in the leading edge of the precursor wave region, the wakefield 
has a sinusoidal/sawtoothed wave form. Since the wakefield in the longitudinal direction is given 
by a static potential, the acceleration and deceleration of electrons repeatedly occurs. Behind 
the sinusoidal/sawtoothed wave train, however, the waveform of the precursor wave is destroyed 
probably due to the nonlinear evolution of the Raman scattering process. The electron is scattered 
by those nonlinear waves and the isotropization of electron are initiated. During this nonlinear 
process, the upstream cold electrons are heated and accelerated. 
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4.1. Particle Acceleration in Leading Edge of Wakefields 

We may simply estimate the maximum attainable energy ffmaxi during this acceleration by 

^maxl ~ eE^akcL-, (15) 

where -Ewake and L are the amplitude of the wakefield and the scale length of the wakefield, respec- 
tively. From Eqs-dHD, ([7]), and p3|) . we have an estimate for the energy of accelerated electrons in 
the downstream frame as 

'"^^^ lALl (16) 



limec^ Y^l + r/a^ 

The efficiency of the electron acceleration increases with increasing upstream bulk Lorentz factor 
71, and the accelerated electron energy is inferred to exceed the upstream bulk energy of ions if 
7i > rrii/me- The energy gain of electrons per particle is sometimes supposed to be limited by the 
ion bulk flow energy per particle, but in this wakefield acceleration due to the radiation pressure, 
it may exceed the ion bulk flow energy in the shock downstream (laboratory) frame. 

We compare the energy of the accelerated electrons between theory and simulation. As we can 
see in Eq. ffTB]) . the energy depends on the amplitude of ao ~ 71 V 2econv {^pe /"^o ) 1 but its dependence 
on (Te is weak. Figure 6 shows the energy spectra at iWpe = 5431 for (ctotai) CTe) = (8 x 10~^,0.4) 
and (3.2 x 10~^, 1.6) with 71 = 10. The format is the same as in Figure 2. Three energy spectra 
were obtained at the same time of ttOpe = 5431. By comparing those energy spectra, one can find 
the maximum energy does not strongly depend on the a parameter. 

Let us qualitatively compare the maximum energy. The solid curve in Figure 7 shows the 
theoretical estimation of the accelerated particle energy in Eq. (jl6p as a function of the upstream 
bulk Lorentz factor 71. The vertical axis is the energy gain normalized by the upstream bulk flow 
energy of electrons. We have assumed cie = 0.1 and econv = 1- This theoretical curve suggests that 
emaxi/7i"^eC^ ~ 8 for 7i = 10, while the simulation result in Figure 1 suggested emax/7i"^eC^ ~ 
20 in the oscillatory wakefield region, and emax/7i"T'eC^ > 50 in the turbulent wakefield region. 
From the energy spectrum in Figure 2, we find that the maximum energy of electrons is about 
emax/7i"ieC^ ~ 270. We think that the theory of ffmaxi = ^E^akeL predicts the particle acceleration 
in the sinusoidal/sawtoothed wakefield region. 

Although there is a small discrepancy between the simulation of Smax/^ifneC^ ~ 20 and the 
theory of £maxi / Jin^eC^ ~ 8, we think that this small difference comes from the non-linear effect 
of the wakefield wavelength shift. As seen in Figure 3, the wavelength of the sawtoothed region is 
longer than that of the tip of the precursor waves. The peak amplitude of the wakefield is shifted 
to k/kp ~ 0.02 — 0.03. By taking into account the change of the wavelength shift with the factor 
of 2 or 3, we can conclude that the agreement between the theory and the simulation is good. 
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4.2. Particle Acceleration in Turbulent Wakefields 

In the turbulent wakefield region, the energy gain is one order of magnitude larger than that in 
the sinusoidal/sawtoothed region. We find that the main energization is occurring in the turbulent 
wakefield region. Let us study the maximum attainable energy in our simulation result. Shown in 
Figure 8 arc the time histories of the maximum energy of electron for Ce = 0.1, 0.4 and 1.6 with 
71 = 10. Their energization occurs in sporadic manner, because the shock upstream does not have 
a time stationary structure but a time variable nature due to the sporadic ejection of the intense 
precursor waves from the shock front. The growth curve is not approximated by a simple linear 
curve, but roughly speaking the acceleration time scale is almost same among three parameters. 
Another more important point is that those curves do not seem to show any saturation behavior 
yet. We have continued the calculation of the time evolution for the case of o"e = 0.1, and we 
find the maximum electron energy rectched at 6jnax 

/7imeC^ ~ 500 at tujpe = 8000 after a couple 
of sporadic energy increase stages. We did not see any signature of the saturation of the growth 
curve. 

The maximum attainable energy in the turbulent wakefield region still remains to be solved 
in the computation study, but we found that the maximum energy is much larger than ffmaxi = 
e-E^wake^- Let us discuss some possible theoretical models to explain such a large energy gain. 



4-2.1. Wakefield Acceleration with Slippage Effect 

In our previous simple model of the particle acceleration of eE'^ake^i we assumed that the 
wakefield structure remains stationary against the particle motion. However, the wakefield is prop- 
agating against the plasma medium, and its phase speed is almost equal to the speed of light. 
Therefore, particles which are being accelerated in the same direction as the wakefield propagation 
can be in resonance with the wakefield electric field. This effect is called as the phase-slippage 
effect. The maximum attainable energy of such a particle in the shock upstream frame can be 
given by, 



nieC^ rrieC^ \c — Vph\ 



(17) 



where the last term of c/|c — Vph\ represents the phase-slippage effect between the particle and the 
wakefield. By using the relationships of the scale length of wakefield L ~ c/upe and the phase speed 
of the wakefield Vph/c ~ (1 — Wpg/wg)^/^ in the shock upstream frame (cf. Appendix B), we obtain, 

2 - 2 ^. (18) 

We have used Cjq 3> Wpe- The slippage effect of the particles with the moving electric field structure 
can be also obtained by the Lorentz transformation of the potential energy. This acceleration 
process is same as that discussed in the laser-beam acceleration (e.g., Tajima and Dawson 1979; 
Esarey et al. 1996). 
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In fact we found particle acceleration due to the slippage effect in our simulation before the 

collapse of the sinusoidal/sawtoothed wakefields, but this does not necessarily mean the production 
of high energy particles in the shock downstream frame, because the accelerated particles have 
negative momenta. The accelerated particles are traveling towards upstream, and cannot directly 
contribute to high energy particles in the shock downstream. To get high energy particles in the 
shock downstream frame, some isotropization/scattering processes in momentum space are needed. 

If we could assume a strong scattering process such that the particles keep the total energy 
but changes the direction of motion in the shock upstream frame, the corresponding energy in the 
shock downstream frame is given by the Lorentz transformation, £:max2 — 7i(l + /3i)£max2 5 and we 
have 

%!-l!elL- (19) 



We have used Coq 271 and oJo/iOpe ~ 0(1). We depicted the theoretical curve by the short 
dashed line in Figure 7 together with the acceleration model of £maxi (solid line) . We have assumed 
that ujo/ujpe = 2.5, rj = 1/2 and econv = 1-0. We think that this optimistic expression probably 
overestimates the acceleration efficiency, because several nonlinear processes that may suppress the 
acceleration efficiency are neglected. The most problematic issue is the scattering of particles, and 
its isotopization process remains an open question. 

Let us think about more realistic scattering process. We think that the isotropization/scattering 
processes occurs during a nonlinear process of the forward Raman scattering (FRS). In the first 
stage of FRS, the injected electromagnetic wave decays into the scattered electromagnetic wave 
and the electrostatic Langmuir wave, and as the second stage of FRS the scattered electromagnetic 
wave can decay again into another electromagnetic wave and the Langmuir wave. In the nonlinear 
FRS, we could expect such a successive Raman scattering process. Since the oscillation frequency 
of the scattered electromagnetic wave (daughter wave) is lower than that of the pump wave, the 
nonlinear Raman can generate the broadband wave spectrum. 

In the upstream plasma frame, the frequency of the initial precursor wave was about u ~ 
271 Wo ~ (4 — 6) X 7iWpe, and due to the above nonlinear FRS the stimulated electromagnetic waves 
are extended down to a; LOpg, namely the broadband spectrum with < a) < (4 — 6) x jiujpe will 
be generated. Those waves satisfy the normal electromagnetic dispersion relation of o)^ = Up^+k^c^. 
The top and left-hand panel in Figure 9 shows schematically the wave dispersion relation of the 
electromagnetic wave in the shock upstream frame, where the plasma is at rest. The thick line 
shows the broadband wave region generated by the nonlinear FRS. 

Similarly, the broadband longitudinal Langmuir wave can be produced in the range of — \/^kp < 
k < —kp and uj = ujpe, where kp = tOpe/c. This electrostatic wave dispersion is shown in the bottom 
and left-hand panel in Figure 9. Note that those electromagnetic and electrostatic waves are 
propagating toward the upstream direction in the upstream frame and uj/k < 0. 

In the shock downstream (simulation) frame, however, a part of electromagnetic and electro- 
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static waves can propagate toward the downstream direction and uj/k > 0. Shown in the right-hand 
panels in Figure 9 is the corresponding wave dispersion relation in the shock downstream frame. By 
virtue of the Lorentz transformation of the wave frequency and the wave number, the dispersion rela- 
tion of Langmuir wave of u) = ujpe in the shock upstream frame is transformed to a; = cjpe/Ti + fcc/3i, 
while the electromagnetic dispersion relation remains the same form of o;^ = LUp^ + k'^cP. The ex- 
isting wave regions are denoted by thick lines, and we can find that the downstream propagation 
of the electromagnetic waves appears for the first quadrant with ujpe < uj < 'jiujpe, and the electro- 
static waves are for the third quadrant with k < and uj < 0. In those broadband spectra in the 
shock downstream (simulation) frame, the largest wave number generated by FRS becomes about 
k ~ jiujpe/c for both the electromagnetic and electrostatic waves. 

Shown Figure 10 is the stack plot of the electromagnetic and the electrostatic wave forms 
around the shock upstream region. The shock front is located around x/{27rc/ujpe) ~ —400. The 
wave forms are superposed from tujpe = 5431 to 5621. In the leading edge of the precursor wave, 
both precursor electromagnetic wave and the electrostatic wakefields are propagating in the negative 
X direction, while near the shock front the downstream propagating waves can be also clearly seen. 
Figure 11 is the wave Fourier spectra in the turbulent wakefield region in the same format as Figure 
3. As approaching toward the shock front, the spectra become broad, but the largest wave number 
was bounded by A; ~ 5 u)pe/c, which is almost equal to (\/3 — l)jiujpe/c for 71 = 10. (Note that 
the Nyquist wave number in our simulation is set to be k^ = tt/Ax = 80.7{ujpe/c).) The results of 
wave propagation behavior in Figures 10 and 11 support the idea of the nonlinear forward Raman 
scattering process. 

Since not only the upstream propagating waves but also the downstream propagating waves 
coexist in the shock upstream region, we can expect that the slippage effect between the particles 
with the positive momenta and the decayed wakefields propagating toward the positive x direction 
in the simulation frame is possible. Let us discuss the maximum attainable energy for this case. 
With the similar discussion in Eq. (jl7p . we get in the shock downstream frame, 



£max3 ^ 6^wake^ C /'9n"\ 
rUeC^ rUeC^ \C-Vph\ 

The acceleration efficiency is again expressed by the product of the wakefield scale L and the 
slippage effect c/|c — Vph\^ but for the case of the nonlinear FRS in the shock downstream frame, 
the scale size L ~ l/Zc increases with decreasing the phase velocity Vph in the range of w/Zc > (see 
the bottom and right-hand panel in Figure 9). We find easily the maximum efficiency appears at 
the largest wave number of ~ (/3i — \/3)7i cjpe/c, i.e., the scale length L ~ c/((\/3 — [3i)^iujpe). 
The slippage effect is approximated by c/|c — Vp}i\ ~ (\/3 — l)^7i. Then the maximum attainable 
energy can be estimated by 

..ax3 ^J^/3-l),a|^^^^,/,^_ ^^^^ 



limec^ y^l + ria^ 

We find that the maximum attainable energy is almost equal to the energy gain in the leading 
edge of the wakefield Emaxi- Although we only discuss the nonlinear process of FRS, but if the 
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backward Raman scattering (BRS) process (see Appendix B) occurs together with FRS, wakefields 

propagating with much higher phase velocities towards the positive x could be generated in the 
shock downstream frame, and then the maximum attainable energy would increase. 

4-2.2. Shock Surfing Acceleration 

Besides the above wakefield acceleration (WFA), wc think that the so-called shock surfing 
acceleration (SSA) can occur in a magnetized shock (e.g. Sagdeev and Shapiro 1973; Lembege and 
Dawson 1984; Katsouleas and Dawson 1985; Ohsawa 1985; Dickman et al. 2000; McClements et 
al. 2001; Hoshino and Shimada 2002; Shapiro et al. 2003). The surfing acceleration is one of the 
direct acceleration mechanisms in plasma, and the particle acceleration occurs during the trapping 
motion due to the electrostatic force. Let us assume the motional electric field Ey in the y direction 
and the ambient magnetic field in the z direction, which mimics the shock upstream region. In 
this situation, the motion of particle is described hy E x B drift, and the plasma is transported 
toward the x direction with the drift speed of cE^/B^. In this situation, for simplicity suppose that 
a stationary electrostatic field structure E^ is induced in a finite zone along the y direction against 
E X B drift motion. The particle can be reflected by the electric force Ex during the gyro-motion, 
and the particle is accelerated along the y direction. However, as long as the Lorcntz force eVyB^/c 
remains smaller than the electric force eE^, the multiple reflection occurs and the particle gains 
energy from the motional electric field Ey. The surfing mechanism is possible when the motional 
electric field Ey exists in the frame of E^ structure at rest. 

In the leading edge of the wakefield, the E^ force balances with the ponderomotive force, 
and the surfing acceleration is not effective. However, behind the sinusoidal/sawtoothed wakefield 
region, the correlation structure between the envelop of the precursor wave amplitude and the 
wakefield is being smeared out, and we could expect the surfing acceleration. As shown in Figure 
1, we find the electrons accelerated toward the y direction in the turbulent wakefield region of 
—600 < x/{2Trc/ujpe) < —400. It should be noted that the signature of SSA in the ion phase space 
at x/{2TTc/LVpe) ~ —580, and the high energy ions seen around x/{27rc/ujpe) ~ —450 are the remnant 
of the SSA accelerated earlier time stage. 

It is known that SSA can provide an unlimited acceleration if Ex > Bz, because the Lorentz 
force {eVyBz < eBz) is always weaker than the electric force (eEx), and because the particle cannot 
escape from the electrostatic potential well (e.g. Katsouleas and Dawson 1985; Hoshino 2001; 
Shapiro et al. 2003). However, in a relalistic shock wave, particles are scattered by non-stationary 
shock behavior and/or turbulent waves, and the unlimited acceleration would not be easily realized. 



-18- 



4-2.3. Combined Acceleration of WFA and SSA 

The particle acceleration in the turbulent wakefield region is the most important in the rel- 
ativistic shock, and we think that not only the wakefield acceleration (WFA) with the slippage 
effect but also the shock surfing acceleration (SSA) is happening. If a successive acceleration could 
occur by switching between WFA and SFA under the turbulent fields, the particle acceleration of 
WFA is not limited by the scale length L described by the nonlinear Raman process, and we might 
adopt the system size Lgyg instead, which is of the order of the scale of an astrophysical object. We 
implicitly assume that the WFA is the main acceleration mechanism and the nonlinear interaction 
of WFA and SSA plays a role for the particle scattering. The maximum attainable energy can be 
estimated by 

eeff , „ eeff ~ -Lsys/(c/Wpe)eeff, (22) 



7imec2 7imec2 ^\ + r/ag 71c 

where represents the efficiency of the combined acceleration. The system size Lgyg can be 
roughly estimated from the propagation speeds of the precursor wave and the shock wave. From 
the Rankine-Hugoniot relation, the shock propagation speed is approximated by c\/r^^, where V 
is the ratio of specific heat. Then the system size is given by 

Lsys - ct(\ - Vr- 1). (23) 

We now compare the above combined acceleration model with the simulation result. From the 
time history of the maximum energy of particle in Figure 8, we may find that emax/(7i'TieC^) ~ 
(0.1 — ^.'2)uj.p(.t. By comparing this with the above theory, we can numerically estimate for the 
acceleration efficiency as 

ee. ~ ^ - \. (24) 

We have assumed F = 3/2, because the plasma thermalization for a perpendicular shock is restricted 
in two-dimensional velocity space perpendicular to the magnetic field. 



5. DISCUSSIONS AND CONCLUSION 

We studied the mechanism of non-thermal particles generation by the energy dissipation of 
the precursor wave in a relativistic shock, and argued that the high energy electrons and ions can 
be accelerated by the electrostatic wakefield induced by the electromagnetic precursor wave by 
an action of the ponderomotive force. We found that the accelerated particle energy can quickly 
attain to £maxi/7i^eC^ ~ 7i in the leading edge of the precursor wave, and we discussed that they 
can be further accelerated during the nonlinear process of the wakefield collapse. We have given 
two theoretical estimations of the maximum energy by taking account of the slippage effect: one 
is £018x2/71^160^ 7i by the strong scattering model and the other is £max3/7i™eC^ 7i by the 
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nonlineax Raman scattering model. The former estimation of £max2 is based on several optimistic 
assumptions such as an isotopization in the upstream plasma frame, and they should be carefully 
discussed further. The latter one of emaxS gave almost the same result as emaxi, but we found 
that the particles can be further accelerated beyond emax3 by the combined nonlinear processes 
of the surfing acceleration and the nonlinear Raman scattering. Then finally we suggested the 
maximum attainable energy as a function of the scale size of an astronomical object Lgys, which 
is £max4/7i"ieC^ ~ Lgys/ {c/u)pe)£eS , where the efficiency of the combined acceleration process is 

£eff ~ 1/6 - 1/3. 

The Wakefield acceleration utilizes the electrostatic field generated in the longitudinal direc- 
tion by the precursor wave, while the conventional shock accelerations such as the diffusive shock 
acceleration, the surfing acceleration and so on uses the motional electric field vx Bi/c < Bi. Since 
the motional electric field is limited to Bi, the maximum energy of the conventional acceleration 
is roughly expressed by eBi x L, where L is the scale size of the system (e.g., Hillas 1984). The 
Wakefield acceleration, however, the amplitude of the longitudinal electric field is not limited by 
the ambient magnetic field Bi, but the amplitude of the electrostatic wakefield -E^ake is of order 
of the electromagnetic precursor wave Bq, which increases with increasing the upstream bulk flow 
energy. Namely, -Ewake ~ Bq. Therefore, the maximum energy of the wakefield acceleration can ex- 
ceed the theoretical limit of the conventional shock acceleration estimated by eBiL. The wakefield 
acceleration becomes important for the case of a low a regime, where the bulk flow energy density 
is larger than the electromagnetic energy density. 

One of the important agents of the wakefield acceleration is plasma composition. If the plasma 
consists of only electron and position, the electrostatic wakefield cannot occur, because the same 
amount of the ponderomotive force works for both positrons and electrons, and because no charge 
separation can be induced. In fact, the large amplitude precursor waves have been observed in 
the electron-positron shock, but the large amplitude i^wakc did not appear in the electron-positron 
shock upstream region by Gallant et al. 1992. In the electron-positron shock, they argued the 
effect of the large amplitude wave into the shock Rankine-Hugoniot relations, and they found that 
the shock structure is slightly modified through an action of the wave pressure of the precursor 
waves for the case of the intermediate cr± ~ 0.1. This is not the wakefield eflFect but just due to the 
ponderomotive force. 

In the multi-component plasma that consists of ion and pair plasma, the wakefield could be 

basically generated. The growth rate of the electrostatic wakefield probably depends on the ratio 
of the ion density to the pair plasma density, and the saturation amplitude of the wakefield may 
be proportional to the density ratio as well. However, due to the self- modulation/feedback process 
between the incoming particles and the wakefield, the wakefield amplitude could be expected to be 
enhanced up to a certain nonlinear saturation level. It is necessary to study carefully the eflFect of 
the plasma composition and the nonlinear saturation amplitude of the wakefield. 

We mainly discussed the electron acceleration of the wakefield, but ions can be also accelerated 
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by the same wakefield. If the electrons can be accelerated to the upstream bulk flow energy of ions, 
7i mjC^, those accelerated electrons become as heavy as the initial upstream ion bulk flow inertia, 
and the acceleration efficiency of the ions becomes same as the accelerated electrons. If the upstream 
bulk Lorentz factor 71 exceeds the mass ratio of ion to electron rrii / rrie , it is expected that both ions 
and electrons are equally accelerated. In this regime, a strong coupling process between electron and 
ion is likely to be important, and not only the electron acceleration but also the ion upstream bulk 
flow modulation may contribute to the shock dynamic evolution (Lyubarsky 2006). The nonlinear 
coupling process between ions and electron should be important and carefully studied. 

The wakefield acceleration due to the radiation pressure of the precursor wave occurs in the 
shock upstream region, and wc expect those pre-accelerated particles could be further accelerated in 
and around the shock front by other shock acceleration mechanisms such as the diffusive Fermi ac- 
celeration and/or the shock surfing acceleration. The wakefield acceleration may serve the diffusive 
Fermi acceleration as the seed particles for further acceleration to much higher energies. 

Although our simulation showed the well-developed structure of the relativistic magnetosonic 
shock, our computational CPU time was not enough to show a fully isotropization process of 
accelerated particles, and a final energy spectrum remains to be solved. On the front side of 
the electromagnetic precursor wave, the oscillatory motion of electrons and ions can be obtained 
in association with the electrostatic wakefield. Behind such oscillatory behavior of plasma, the 
wave form of the wakefield quickly collapses into a turbulent state, and electrons and ions start 
to become isotropized. During this isotropization, particles can be further accelerated. We may 
expect the stochastic particle acceleration with a random accelerating-decelerating phase in the 
nonlinear collapse region of the wakefields, which might result in the formation of a power-law energy 
spectrum (e.g., Chen et al., 2002). The nonlinear processes that control a long time relaxation of 
the accelerated particles remain an open issue. 

Our simulation study is only a first step towards the understanding of the rich potential 
wakefield acceleration in a relativistic shock. It is important to study further our simple acceleration 
process from several other aspects: one is the radiation loss processes such as the synchrotron loss 
and Compton scattering. Another is the wave coherency of the precursor wave that is probably 
the most important unresolved subject, because the coherency is required for the ponderomotive 
force (cf. Appendix A). Otherwise the effect of the radiation pressure force enters into the standard 
Thomson scattering regime, which interaction force is much weaker than the ponderomotive force 
in collective plasma phenomena. 

I am grateful to T. Amano, C. Jaroschek, Y. Lyubarsky, M. Lyutikov, K. Mima , K. Nagata, 
T. Tajima and Y. Takagi for useful discussions. 
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A. PONDEROMOTIVE FORCE BY PRECURSOR WAVES 

Let us quickly review the ponderomotive force in plasma. Consider a spatially inliomogeneous 
electromagnetic wave at frequency w, 

E{x,t) = Eo{x,t)cos{u;t), (Al) 
B{x,t) = Bo{x,t)sin{ut), (A2) 

where E{x,t) and B{x,t) satisfy the Faraday's law, 

B{x,t) = -c [ V X Eo{x,t)cos{ujt)dt. (A3) 
Jo 

We have assumed that the variation of the wave amplitude in x direction, and the transverse, 
electromagnetic wave in the y — z plane. In addition to the coherent wave with the frequency w, 
we have also assumed that the modulation of the wave can be described by EQ{x,t) and BQ(x,t). 

The equation of motion of a charged particle in this field can be written 

dv 1 -* V -* \ 

m— = e Eo{x, t) cos{ut) + - x B{x, t) , (A4) 
at \ c J 

It is expected that the motion of the charged particle has two time scales: one is a simple harmonic 
motion in response to a high frequency, oscillating electric field, and the other is a slow movement 
of the center of oscillation. To understand this, we introduce the variable of v as 

^ = U + u (A5) 

where U is slowly varying quantity, while u is rapidly varying quantity with an oscillating electro- 
magnetic field. We have assumed that \U\ <S l^l and cEq/ {muic) <^ 1. To the first order, we have 
the fast oscillating motion of 

du 6 -* 

— = —Eo{x,t)cos{Lot) (A6) 
By taking into account the constrain of (u) = 0, the above solution can be given by 

e /•* - 

■=— / Eo(x,t) cos( iot)dt (A7) 
m Jo 

To the next order, we will seek the slowing varying motion of the charge particle. The oscillation 
average yields, 

dU _ e -i. , , 

m— = -{ux B) (A8) 

= ( / Eo{x,ti) cos{u;ti)dti X / V x Eo{x,t2) cos{ut2)dt2) (A9) 

nT- Jo Jo 

If we could assume Eo{x,t) = Eo{x) and Bo{x,t) = Bo{x), we obtain 

dU 



u 
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Namely 

dU _ , / * -, , X 

m— = -eVcppond, (All) 

where 

'/'pond = J^£;o'(a;)- (A12) 

We find that particle can be accelerated due to the wave pressure force, which is called as the "pon- 
deromotive force" . The ponderomotive force is proportional to the gradient of the wave pressure, 
and is independent of the sign of the charge, but the force of electron is much larger than that of 
ion. 

In the above discussion we assumed "coherency" of the electromagnetic wave, i.e., Eo(-J',t) = 
Eq(x) and BQ{x,t) = Bq{x). However, if there is no correlation between EQ{x,ti) and EQ(x,t2) for 
two different times of ti and t2, the above discussions is no longer valid. Let us assume that two 
point correlation of the electromagnetic field can be described as 

{Eo{x,h) ■ Eo{x,t2)) = {El{x,t)) exp{-K^{h - t2f). (A13) 

Similarly we may also assume that the correlation between the irregular part of E(){x,t) and the 
coherent part of cos{ut) is destroyed within the time scale of k, i.e., 

Eo{x,t)cos{ujt)dt ~ (£;^(x,0>^^^e-'^*^^^^. (A14) 
Jo ^ 

In this case the ponderomotive force can be modified as 

We can find that the ponderomotive force is effective as long as > k. For a high frequency 
electromagnetic (light) wave with uj ^ ujpf., the correlation of the wave phase could be easily 
satisfied in most of collective plasma processes. 



B. RAMAN SCATTERING IN INTENSE ELECTROMAGNETIC WAVES 

Under a nonlinear plasma process, it is known that a large-amplitude, electromagnetic light 
wave decays into a relativistic Langmuir wave plus a scattered electromagnetic light wave by the 
stimulated Raman scattering instability (e.g., Mima and Nishikawa 1988; Kruer, 1988). When the 

frequency and wave number of the pump electromagnetic light wave (photon) is given by (wq, ko), 
the scattered electromagnetic light wave (photon) with (wi , ki ) and the induced Langmuir wave 
(plasmon) with {L02,k2) have to satisfy the following matching condition. 



CJQ = LOI+UJ2, 

ko = ki + k2, 



(Bl) 
(B2) 
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and the strong coupling can occur when those waves satisfy the normal modes of the responsive 
plasma dispersion relation. Namely, the electromagnetic waves satisfy the dispersion relation of 
ujf = + kfc^ where i = or 1, and the Langmuir wave is = oj^^. For simplicity, we neglected 
effects of a finite plasma temperature and magnetic field. Also we assumed the reference frame 
where the ambient plasma is at rest in this Appendix. 

The Raman scattering instability can be understood as follows: let us assume that an electro- 
magnetic (light) wave with the electric field Et is propagating through a plasma medium with the 
density fluctuation 5n associated with a Langmuir wave. Since the electrons are oscillating in the 
electromagnetic (light) wave with the velocity vt-, the transverse electric current SJt = —evTSn 
can be generated. If the nonlinear coupling of vt and 6n produces the proper frequency and wave 
number which matches with the normal mode of the electromagnetic (light) wave, this scattered 
electromagnetic (light) wave with the electric field 6E interferes with the incident light wave to 
produce a variation in the wave pressure, namely, VE"^ = 2V{ElSE). If we assume the frequency 
of the incident electromagnetic (light) wave u>o and the frequency of the density fluctuation L02, 
the induced electric current SJ and the electric field SE may have the frequency of ui = loq — lj2- 
Then the frequency of the ponderomotive force has luq — loi, which is same as the original density 
perturbation of the frequency with uj2- Due to this feed-back process through the ponderomotive 
force, the plasma density can be enhanced. This is an essence of the stimulated Raman scattering 
instability. 

There are basically two possible Raman scattering processes in plasmas: one is the forward 
Raman scattering (FRS) in which a scattered electromagnetic wave propagates in the same direction 
as the pump electromagnetic wave. The other process is the backward Raman scattering (BRS) in 
which the propagation of a scattered electromagnetic wave is opposite to the pump electromagnetic 
wave. 

For the case of FRS, if the pump wave frequency loq is much larger than the plasma frequency 
Upe, we can easily show that the wave length of the Langmuir wave k2 becomes close to the electron 
inertial length, 

^2 ~ — , (B3) 

and the phase velocity of the Langmuir wave (plasmon) 102/ k2 can also be estimated from the 
matching condition of Eqs.(Bl) and (B2), 

iJ2 ^ ^ _ ^ _ 

k2 Y ooq dko 

We find that the phase velocity of the Langmuir wave approaches the group velocity of the pump 
electromagnetic wave, which is almost equal to the speed of light c. 

For the case of BRS, the Langmuir wave in general has a larger wave number than the pump 
electromagnetic wave. If the electron plasma temperature is finite, the Langmuir waves (plasmon) 



(B4) 
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with a large wave number may suffer from Landau damping. Therefore, BRS may play a minor 
role on the nonlinear wave process of the Raman scattering (e.g., Kruer 1988). 

If we take into account a uniform magnetic field perpendicular to the wave vector in the shock 
precursor region, the electromagnetic waves should read the extraordinary mode (e.g., Hoshino and 
Arons 1991), and its dispersion relation is given by 





(B5) 



For the case of a small a the effect of magnetization may be negligible. Furthermore, the dis- 
persion relation of the Raman scattering instability has essentially the same form as that for an 
unmagnetized plasma (e.g., Mima and Nishikawa 1988). 



-25- 



REFERENCES 

Amano, T., & Hoshino, M., 2007, ApJ, 661, 190 
Asseo E. and Sol, H. 1987, Phys. Rep., 148, 301 

Armstrong, T. P., Pesses, M. E., & Decker, R. B., 1985, in Collisionlcss Shocks in the Hcliosphere: 
Review of Current Research, eds. B. T. Tsurutani & R. G. Stone (Washington, DC, American 
Geophysical Union), 271 

Axford, W. L, Leer, E., & Skadron, G. 1977, Proc. 15th Internat. Cosmic Ray Conf., 11, 132 

Bauer, D., Mulser, P., & Steeb, W.-H. 1995, Phys. Rev. Lett., 75, 4622 

Begelman, M.C., Blandford, R.D. and Rees, M.J., 1984; Rev. Modern Phys., 56, 255 

Bell, A. R. 1978, MNRAS, 182, 147 

Bieber, J.W., Matthaeus, W.H., Smith, C.W., Wanner, W., Kallenrode, M.B., and Wibberenz, G., 
1994, ApJ, 420, 294 

Blandford, R. E. & Ostriker, J. P. 1978, ApJ, 221, L29 

Blandford, R., & Eichler, D. 1987, Physics Report 154, 1, 1 

Carilli, C. L., & Barthel, P. D. 1996, Astron Astrophys. Rev., 7, 1 

Cesarsky, C. J., 1992, Nucl. Phys. B (Proc. Suppl.), 28, 26 

Chen, P., Tajima, T., & Takahasi, Y. 2002, Phys. Rev. Lett., 90, 161101 

Coroniti, F. V. 1990, ApJ, 349, 538 

Dieckmann, M. E., McClements, K. G., Chapman, S. C, Dendy, R. O., & Drury, L. O'C. 2000, 
Astron. Astrophys., 356, 3771 

Drury, L O'C. 1983, Rep. Prog. Phys., 46, 973 

Esarey, E., Sprangle, P., Krall, J., & Ting, A. 1996, IEEE Trans. Plasma Sci., 24, 252 
Gallant, Y. A., Hoshino, M., Langdon, A. B., Arons, J., & Max, C. E. 1992, ApJ, 391, 73 
Hillas, A. M., 1984, Ann. Rev. Astron. Astrophys. 22, 425 
Hess, V. 1912, Physik. Zeitschr., 13, 1084 

Hededal, C. B., Haugboole, T., Frederiksen, J.T., and Nordlund, A. 2004, ApJ, 617, L107 
Hoshino, M. & Arons, J. 1991, Phys. Fluids, B3, 818 



-26- 



Hoshino, M., Arons, J., Gallant, Y. A., & Langdon, A. B. 1992, ApJ, 390, 454 
Hoshino, M. 2001, Prog. Theor. Phys. Suppl., 143, 149 
Hoshino, M. & Shimada, N. 2002, ApJ, 572, 880 

Jaroschek, C. H., Treumann, R. A., Lesch, H., & Scholer, M. 2004, Phys. Plasmas, 11, 1151 
Jones, F.C. and Ellison, D.C., 1991, Space Sci. Rev., 58, 259 
Junor, W., Biretta, J. A., & Livio, M. 1999, Nature, 401, 891 
Kato, T. 2005, Phys. Plasmas, 12, 080705 

Katsouleas T., & Dawson, J. M. 1983, Phys. Rev. Lett., 51, 392 
Kirk, J. G., & Skjaeraasen, P. 2003, ApJ, 591, 366 

Koyama, K., Petre, R., Gotthelf, E. V., Hwang, U., Matsuura, M., Ozaki, M., & Holt, S. S. 1995, 
Nature, 378, 225 

Kruer, W. L., 1988, The Physics of Laser Plasma Interaction, Frontiers in Physics, Addison- Wesley 
Publ. 

Langdon, A. B., Arons, J. & Max, C. E. 1988, Phys. Rev. Lett., 61, 779 
Lembege, B. & Dawson, J.M. 1984, Phys. Rev. Lett., 53, 1053 
Lyubarsky, Y. E., 2006, ApJ, 652, 1297 

McClements, K. G., M. E. Dieckmann, Ynnerman, A., Chapman, S. C., k. Dendy, R. O. 2001, 
Phys. Rev. Lett., 87, 255002 

Meisenheimer, K. Roser, H.-J., Hiltner, P. R., Yate, M. G., Logair, M. S., Chini, R., & Perley, R. 
A. 1989, A&A, 219, 63 

Mima, K. &: Nishikawa, K. 1988, Parametric Instabilities and Wave Dissipation in Plasmas, in 
Basic Plasma Physics, eds. M.N. Rosenbluth and R.Z. Sagdeev, North-Holland Phys. Publ. 

Mourou, G.A., Tajima, T. and Bulanov, S.V., 2006, Rev. Mod. Physics, 78, 309 

Nishikawa, K. -I., Hardee, P., Richardson, G., Preece, R., Sol, H., k, Fishman, G. J. 2003, ApJ, 
595, 555 

Ohsawa, Y. 1985, Phys. Fluids, 28, 487 

Robinson, P. A., 1997, Rev. Modern Physics, 69, 507 

Romanova, M. M. k Lovelace, R. V. E. 1992, A&A, 262, 26 



-27- 



Sagdeev, R. Z. & Shapiro, V. D. 1973, JETP Lett. Engl. Transl., 17, 279 

Shapiro, V.D. and Ucer, D., 2003, Planet. Space Sci., 51, 665 

Sonnerup, B. U. O., 1969, J. Geophys. Res., 74, 1301 

Sprangle, P., Esarey, E. & Ting, A. 1990, Phys. Rev. Lett., 64, 2011 

Tajima, T., & Dawson, J.M. 1979, Phys. Rev. Lett., 23, 267 

Thompson, C. 1994, MNRAS, 253, 141 

Thompson, C. & Duncan, R. C, 1995, MNRAS, 275, 255 

Vietri, M. 1995, ApJ, 453, 883 

Yang, T.Y.B., Arons, J., and Langdon, A.B. 2004, Phys. Plasmas, 1, 3059 

Zenitani, S. & Hoshino, M. 2001, ApJ, 562, L63 

Zenitani, S. & Hoshino, M. 2005, Phys. Rev. Lett., 95, 095001 



This preprint was prepared with the AAS IM^^X macros v5.2. 



- 28 - 




SDace olots — X and — X. The second Danel from the bottom is the electrostatic fields E.,. 



-29- 




■I 4; t f 



T — r — I — T — I — I — r—i — r — p"^T — i — r — i — i — i — r-n — i — | — r^-i — i — j-n — i — r — r— i — | — i — i — i — i — r-T — t — i — r 




IQO 200 igft 4Qp, 

■rhi 



Fig. 2. — Energy Spectra for ions and electrons in the shock upstream (dashed hnes) and the shock 
downstream (sohd hnes). 




Fig. 3.— 



Waveform and wave spectra for and in shock upstream region. 
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Fig. 4. — Precursor wave energy density as a function of a±. (a) Wave energy normalized by the 
upstream ambient magnetic energy, and (b) wave energy normalized by the upstream bulk flow 
energy. 
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Fig. 5. — Comparison of the upstream wave amplitudes between the wakefield acceleration theory 
and the PIC simulation. The horizontal axes are ae, while the vertical axes are the normalized 
amplitudes: (a) the electromagnetic precursor wave Bq, and (b) the electrostatic wakefield -E^ake- 




Fig. 6. — Ion and electron energy spectra for cje = 0.4 (left) and = 1-6 (right) in the same format 
as Figure [2j 
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Fig. 7. — Two different models of the maximum attainable energy of accelerated particles. The 
solid line (emaxi) stands for a simple acceleration model in the leading edge of the wakefields, and 
the dashed line (£max2) is the acceleration model with the slippage effect and the strong particle 
scattering. 




Fig. 8. — Time history of maximum energy. The sohd hne, long dashed hne, and short dashed hne 
stand for cTg = 0.1, 0.4 and 1.6, respectively. 
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Fig. 9. — The dispersion relation in — A; space for both electromagnetic precursor waves (top) and 

electrostatic wakcfields (bottom). The generated waves due to the nonlinear Raman scattering are 
shown by thick lines. The shock upstream frame (left) and the downstream frames (right). 
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Fig. 10. — Stack plot of electrostatic (top) and electromagnetic (bottom) wave forms for 
o"e = 0.1 and tOpeto = 5431. 
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Fig. 11. — Waveform and wave spectra for and in the turbulent wakefield region. Same 
format as Figure 3. 



